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We present an introduction to the principles behind atomic energy level calculations with Quan- 
tum Electrodynamics (QED) and the two-time Green's function method; this method allows one 
to calculate an effective Hamiltonian that contains all QED effects and that can be used to predict 
QED Lamb shifts of degenerate, quasidegenerate and isolated atomic levels. 



INTRODUCTION 



This contribution is concerned with the evaluation of atomic energy levels with QED. Such an evaluation yields 
stringent tests of QED in strong electric fields, whereas (/-factor experiments and calculations currently probe QED 
in situations where the magnetic field can be treated perturbatively 

The nuclear Coulomb field experienced by the inner levels of highly- charged ions makes the electrons reach relativistic 
velocities. Such simple physical systems are thus particularly interesting for testing relativistic effects in quantum 
systems (for example, see Refs. [[j]]^] for experimental results with lithiumlike ions). Theoretical predictions of energy 
levels in such systems obviously require the use of QED. 

Experiments have reached an accuracy that shows that extremely accurate evaluations of QED effects are also 
needed in helium. Experiments performed during the last ten years in the spectroscopy of this atom have become two 
orders of magnitude more precise than the current theoretical calculations (see for instance Refs. ||^|] and references 
therein). 

Several experiments are now focusing on helium and heliumlike ions, and especially their ls2p 3 Pj fine structure |^-(^] 
such experiments have implications in metrology, as they could provide a measurement the fine structure constant and 
provide checks of theoretical higher-order effects. Very precise theoretical calculations of energy levels in heliumlikc 
ions can be also important in the investigation of parity violation || . 

Predictions of energy levels are usually more difficult to obtain for states with one or more open shells (retardation 
in the interaction and exchange of electrons must be included, and there can be quasidegenerate levels). Only a few 
calculations of excited energy levels in heliumlike and lithiumlike ions have been performed up to now; the first results 
have been published quite recently |Tc|-|l^ ]. In regards to QED shifts of quasidegenerate levels, they have only been 
obtained this year for the first time |fl| , with the help of the method that we present in this talk. 

THEORETICAL METHODS 

As is well known, relativistic electrons orbiting a nucleus are well treated with the Dirac equation, in which the 
nucleus can be considered as point-like or not. We thus treat the binding to the nucleus non-perturbatively by using 
"Bound-State QED" jlj^,[l4j (the coupling constant of the nucleus-electron interaction is Za, which is not small for 
highly-charged ions). In this formalism, however, QED effects are taken into account by treating the electron-electron 
interaction perturbatively (with coupling constant a), and both the electron and photon fields are quantum fields (i.e., 
in second-quantized form); the only difference with the free-field case used in high-energy physics is that electronic 
creation and annihilation operators create and destroy atomic states instead of free particles. 

A few methods allow one to extract energy levels from the Bound-State QED Hamiltonian: the two-time Green's 
function method Jl5|-p^|, the method being developed by Lindgren (based on Relativistic Many-Body Perturbation 
Theory merged with QED) fl8|,|l9|| , the adiabatic S'-matrix formalism of Gell-Mann, Low and Sucher [^o|, and the 
evolution operator method J2l| , [22| ]. Some other methods yield atomic energy levels, but they include QED effects 
only partly or approximately (such as the multiconfiguration Dirac- Fock method J23|| , configuration interaction calcu- 
lations Q and relativistic many-body perturbation theory |25|]), 

However, only two methods can in principle be employed in order to calculate energy levels of quasidegenerate 
atomic states [e.g., the (2s2p!/ 2 )i and the (2s2p 3 / 2 )i levels in heliumlike ions, which are experimentally important]: 
the two-time Green's function method and the method being elaborated by Lindgren. We present in this talk a non- 
technical introduction to the first method. The two-time Green's function method has also the advantage of yielding 
a simpler renormalization procedure than the Gell-Mann-Low-Sucher method in the case of degenerate levels pG.E7| . 
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THE TWO-TIME GREEN'S FUNCTION METHOD 



All the methods that extract atomic energy levels from the Bound-State QED Hamiltonian study the propagation 
of electrons between two different times. The methods differ in the number of infinite times used: 

(a) in the Gell-Mann-Low-Sucher method, the atomic state under consideration evolves from time — oo to time 
+00 with an adiabatic switching of the interaction; (b) in Lindgrcn's formalism ]l8| , |l9| |. the evolution is from time 
— 00 to time 0, which avoids problems associated with the two infinite times in the S'-matrix approach of Gell-Mann- 
Low-Sucher; (c) in the two-time Green's function method, that we present here, the adiabatic switching is completely 
avoided by studying the propagation of electrons between two finite times. We note that adiabatic switching of the 
interactions is physically motivated in the study of collisions between particles that start very far from each other, 
but this switching is not so easily related to the physical description of the orbiting electrons of an atom. 



The Green's function 




Time f Time f 



FIG. 1. The two-time Green's function represents the probability amplitude for going from one position of the electrons to 
another position. 

The effective Hamiltonian derived from QED by the two-time Green's function method has matrix elements be- 
tween the various degenerate and/or quasidegenerate states under study; the eigenvalues of this Hamiltonian are the 
atomic energy levels predicted by QED (to a given order). This effective Hamiltonian is however not associated to a 
Schrodinger equation of motion; our Hamiltonian is equivalent to the submatrix used in the perturbation theory of 
degenerate and quasidegenerate states; in this respect, the approach of the two-time Green's function method differs 
from the spirit of the Bethe-Salpeter equation. 

The QED Hamiltonian of the method is defined with the help of a Green's function that represents the propagation 
of N electrons between two different (finite) times (N is the number of electrons of the atom or ion that we want to 
study); this propagation is represented in Fig. |l|. 



Atomic energies are in the Green's function 



|gN=2 (E) | >, 
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Poles at the energy 
levels of a 2-electron ion 
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FIG. 2. The two-particle Green's function as a function of energy contains all the information about the atomic energy levels 
of a two-electron atom or ion. 
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The energy levels of an ./V-electron ion or atom can be recovered by studying the energy representation Q N (E) 
of the Green's function, i.e., by doing a Fourier transform: this function has (simple) poles at the atomic energy 
levels ]l5|— |lTt] . Such a result is similar to the Kallen-Lehmann representation p8|| . As an example, Fig. || depicts the 
poles of the two-particle Green's function. 

The two-time Green's function method provides a way of mathematically extracting from the Green's function the 
positions of the poles, i.e., the atomic energy levels Jl^]; the procedure handles degenerate and quasidegenerate atomic 
levels without any special difficulty |29| . One of the basic ideas behind the pole extraction is found in the following 
mathematical device, which uses any contour Tq that encloses the pole in order to find its exact position: if the 
function g(E) has a simple pole at E — Eq, then we have from complex analysis 

& dEE x g(E) 

E = J -^j ; (1) 

f dEg(E) 

the contour Tq is only required to encircle the pole and to be positively oriented, as shown in Fig. ||. Since the Green's 
function has simple poles at the atomic energy levels jlj] , Eq. (]l|) is a way of obtaining them. 




FIG. 3. The exact atomic energies can be recovered through a contour integration of the Green's function. 

When QED shifts of degenerate of quasidegenerate levels are calculated, the scalar Green's function g of Eq. ([j]) is 
simply replaced by a finite-size matrix defined on the space of levels under consideration [EdI . 



GRAPHICAL CALCULATIONS 



Obviously, analytic properties of the Green's function p7| are important in the evaluation of Eq. (|l|). We have 
developed a set of graphical techniques that allow one to obtain the Laurent series of the Green's function Q N (E) by 
a systematic procedure. The idea behind these techniques consists in displaying the analytic structure of the Green's 
function step by step; each step explicitly extracts one singularity, and we proceed until we have exhausted all the 
singularities of the Green's function; at this point, contour integrals such as Eq. (Q) can be calculated quite simply. 

It is impossible to give here a full account of the method we use for deriving the effective, finite-size QED Hamilto- 
nian. However, we can mention a particular feature of our calculational strategy: a very special "particle" appears in 
our algorithm; this particle is quite simple since it "disintegrates" immediately (zero life time) and cannot move (zero 
probability for going from one position to a different one). In mathematical terms, the coordinate-space propagator 
of this particle is a four-dimensional Delta function 5^ [(x, t); (x',t')] that we represent by a special line in Feynman 
diagrams. 



THE SCREENED SELF-ENERGY 



The experimental accuracy on transition energies is so high that second-order (i.e., two-photon) effects must be taken 
into account in order to compare experiments with theory. We thus have very recently calculated the contribution of 
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the self-energy screening []30| to the QED effective hamiltonian; this contribution corresponds to the following physical 
processes: 




Our result is part of the current theoretical effort developed with the aim of matching experimental precisions. 
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